/*******************************************************************************
*   			Impacts of a Large-Scale Parenting Program: 				   *
*					Experimental Evidence from Chile						   *
********************************************************************************


	REQUIRES:	"${dt_final}panel_RR_indexes.dta"
	CREATES:	${out_tables}/T16A_selec_unobs.xls
				
	WRITEN BY:  Italo Lopez  [italolop@usc.edu]

********************************************************************************
	Prepare data
*******************************************************************************/

	* Load data
	use "${dt_final}/panel_RR_indexes.dta", clear


	* Establish macros for regressions

/*******************************************************************************
*   Set macros for regressions
*******************************************************************************/
	
	* Outcomes
	*-----------------------------------------------------------
	local childoutcomes  "zDCCSPunt zptevir_irtscore2 zt_score zcbcl_t_e zcbcl_t_i zibatt_sspi zibatt_ssai zibatt_sssr zas_raw_score" 
	local behaviors "home_score se_score disc_score" 
	local beliefs "zpstyle1 zpstyle2 zpstyle3 zPSCS zPSSS_family zPSSS_friends zPSSS_others zPSI_Distress zCESD zpacotis"	
	local child "zptevir_irtscore2 zcog_index zsoc_index" 

	* Controls
	*-----------------------------------------------------------

	macro  def Xvar1 "i.age_year i.gender"  	
	macro  def Xvar2 "i.age_year i.gender inc_qtaut_old i.pc_edu_mdsfin hh_mem hh_tipo"  	
	macro  def Xvar2old "i.gender inc_qtaut_old i.pc_edu_mdsfin_old hh_mem hh_tipo"  

	
	* Parental characteristics: WAIS and BFI and imputed wais 
	*-----------------------------------------------------------

	macro def parent_cha4 "BFI_ext BFI_ope BFI_con BFI_neu BFI_agr wais2 flag_wais"
	

/*******************************************************************************
	Table A16: Selection on Unobservables 
*******************************************************************************/


local intermediate "zfbeh_gsem zfbel_gsem"

foreach var of local intermediate {				
				*treated nep B
				xi: xtreg `var' i.TIPO  $Xvar1 if  treat2==1 & x==1, fe i(CENTRO_SALUD)
				estimates store `var'_b
				test _ITIPO_2 _ITIPO_3
				estadd scalar p1=r(p)
				
				*treated nep I
				xi: xtreg `var' i.TIPO  $Xvar1 if  treat3==1 & x==1, fe i(CENTRO_SALUD)
				estimates store `var'_i
				test _ITIPO_2 _ITIPO_3
				estadd scalar p1=r(p)
				
				*untreated
				xi: xtreg `var' i.TIPO  $Xvar1 if  treat2==0 & treat3==0 & x==1,  fe i(CENTRO_SALUD)
				estimates store `var'_u	
				test _ITIPO_2 _ITIPO_3
				estadd scalar p1=r(p)
		
				
				
			}
			

			#delimit ;
			
		
esttab zfbeh_gsem_b zfbeh_gsem_i zfbeh_gsem_u  using "${out_tables}/appendix/T16A_selec_unobs.csv", replace f 
	 keep(_ITIPO* ) 
	 cells("b(fmt(3)star)" "se(fmt(3)par)") 
	 star(* 0.10 ** 0.05 *** 0.01) 
	 varlabels(n_ses_b "NEP-B" n_ses_i "NEP-I") 
	 title("Panel A: Parental Behavioral Index")  
	 stats(p1 N, fmt(3 0) 										  
	 labels("P-value NEP-B=NEP-I" "Observations" )) legend;
	
	
esttab  zfbel_gsem_b zfbel_gsem_i zfbel_gsem_u using "${out_tables}/appendix/T16A_selec_unobs.csv", append f 
 	keep(_ITIPO* )
 	cells("b(fmt(3)star)" "se(fmt(3)par)") 
 	star(* 0.10 ** 0.05 *** 0.01) title("Panel B: Parental Beliefs Index") 
	varlabels(n_ses_b "NEP-B" n_ses_i "NEP-I") 
	stats(p1 N, fmt(3 3 0) 
	labels("P-value NEP-B=NEP-I" "Observations")) legend;
	
